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Abstract 



Several different methods have recently been proposed for calculating the 
motion of a point particle coupled to a linearized gravitational field on a curved 
background. These proposals are motivated by the hope that the point par- 
ticle system will accurately model certain astrophysical systems which are 
promising candidates for observation by the new generation of gravitational 
wave detectors. Because of its mathematical simplicity, the analogous system 
consisting of a point particle coupled to a scalar field provides a useful context 
in which to investigate these proposed methods. In this paper, we generalize 
the axiomatic approach of Quinn and Wald in order to produce a general 
expression for the self force on a point particle coupled to a scalar field follow- 
ing an arbitrary trajectory on a curved background. Our equation includes 
the leading order effects of the particle's own fields, commonly referred to as 
"self force" or "radiation reaction" effects. We then explore the equations of 
motion which follow from this expression in the absence of non-scalar forces. 

I. INTRODUCTION 



There has been much recent interest in calculating the motion of astrophysical systems 
which emit gravitational waves in anticipation of data from a new generation of detectors. 
Full three-dimensional numerical simulations are required in order to produce useful results 
for many of the most promising observational candidates, such as colliding black holes. 
However, there also exists a large class of systems which can be accurately modelled by a 
small isolated body moving in the fixed background created by a much larger body (e.g., a 
solar mass star falling into a supermassive black hole). For such a system, we might hope to 
produce useful results by treating the smaller object as a point particle and introducing the 
effects of its fields and internal structure as perturbations to the background geodesic orbit. 

The perturbations due to the particle's own fields, commonly called "radiation reaction" 
or "self force" effects, are particularly important because they include the forces responsible 
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for the decay of the body's orbit. If both the background spacetime and the unperturbed 
orbit of the body possess enough symmetry, it is possible to infer the effects of these forces 
on the orbit from global conservation principles: one calculates the energy and/or angular 
momentum radiated to infinity by a particle in geodesic motion, and then modifies the orbit 
to reflect this energy and angular momentum loss in a time-averaged fashion. (Obviously, 
this procedure can be iterated if greater accuracy is required.) Some justification for this 
method is provided by Quinn and Wald |T|. However, in the absence of such symmetries, 
it is necessary to directly calculate the effects of the local fields in the neighborhood of the 
particle. Unfortunately, this problem is ill-posed, since the fields diverge in the neighbor- 
hood of the particle's world line, so that any such local calculation must include a rule for 
extracting the appropriate finite part of these divergent fields. 

There is an extensive literature devoted to this regularization problem. In 1938, Dirac |2| 
reproduced the force expression (originally given by Abraham ||) for a point particle coupled 
to an electromagnetic field in Minkowski spacetime by imposing local energy conservation on 
a tube surrounding the particle's world line and subtracting the infinite contributions to the 
force through a "mass renormalization" scheme. In 1960, Dewitt and Brehme generalized 
this approach to an arbitrary curved background spacetime. (A trivial calculational error in 
their paper was later corrected by Hobbs ||.) More recently, Mino et al. || further adapted 
this approach to produce a force expression for a point particle coupled to a linearized 
gravitational field on a vacuum background spacetime, and Quinn and Wald |7| rederived 
both the electromagnetic and gravitational forces using an axiomatic approach which, in 
effect, regularizes the forces by comparing forces in different spacetimes. 

There has emerged from this work a consensus regarding the correct equation of motion 
for a particle coupled to electromagnetic fields on an arbitrary curved background and for a 
particle coupled to linearized gravitational fields on a vacuum background. In principle, the 
latter equation allows one to calculate the dynamics of the astrophysical systems of interest 
described above. In practice, however, very little progress has been made in applying either 
equation of motion to concrete physical examples for two reasons. First, given a world line 
in an arbitrary spacetime, the calculation of the associated retarded fields is a complex and 
difficult problem. Second, once these fields are calculated, the equations of motion require 
one to identify that portion of the retarded field at each point of the world line which arises 
from source contributions interior to the light cone. This part of the field is often called the 
"tail term," and most approximation schemes for calculating the retarded field entangle the 
tail and non-tail contributions to the field. 

Nevertheless, some progress has been made, notably in the electromagnetic case. In 
1964, DeWitt and DeWitt || calculated the tail term for an electromagnetic particle in a 
circular orbit on a Schwarzschild background to leading order in the background curvature 
and the velocity of the particle. In 1980, Smith and Will H calculated the force on an 
electromagnetic particle held static on a Schwarzschild background, essentially by repeating 
DeWitt and Brehme's local stress-energy conservation argument. Neither result has been 
generalized to the case of a massive particle coupled to gravitational fields, nor has there been 
any direct progress on the more complex systems which are of interest to the gravitational 
wave astronomy community. However, several new ideas have emerged in recent years which 
may lead to further progress. Ori |10[ has suggested an alternative regularization scheme 
involving averaging of multipole moments which is better adapted to concrete calculations, 
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while others have suggested a hybrid scheme in which the tail term is calculated through a 
combination of Hadamard expansion techniques for small distances and multipole techniques 
for larger distances [O . 



It is clearly important to test these ideas. In particular, we must know whether these 
schemes are equivalent to the equations of motion discussed above. Because of its mathe- 
matical simplicity, one natural system in which to explore all of these questions is that of a 
point particle coupled to a scalar field. Motivated by this, several researchers have begun to 
apply the ideas discussed above to the scalar system. In particular, Ori's method has been 
applied to the motion of scalar particles in the Kerr spacetime |I3 and in the Schwarzschild 



spacetime [p~3| , |T4|] , and Wiseman ||T5[ has adapted the calculation of Smith and Will in order 
to calculate the force on a scalar particle held static in the Schwarzschild spacetime. In 
the present paper, we generalize the axiomatic approach of Quinn and Wald M in order 
to produce the general equation of motion for a point particle coupled to a scalar field on 
an arbitrary background spacetime. It is hoped that this general expression will be useful 
in evaluating the validity of the calculational schemes described above for the scalar case, 
and that this comparison will ultimately help to clarify the relationship between the various 
methods which have been proposed for the electromagnetic and gravitational cases. 

In Sec. |TJ, we derive an expression for the force on a particle following an arbitrary 
trajectory in curved spacetime. Then, in Sec. |TTT| , we explore the equations of motion which 
follow from this expression in the absence of non-scalar forces. 



II. THE SCALAR FORCE 

Given a spacetime containing a particle world line and a Klein-Gordon field sourced by 
the particle, we wish to define the total scalar force f§ on the particle at each point of the 
world line, including so-called self-force or radiation reaction effects. For an electromagnetic 
point particle in flat spacetime, an expression of this sort was first given by Abraham [|J in 
1905, was later rederived in a relativistic context by Dirac 0, and is often found in textbooks 
(e.g., Jackson [0). However, since there are no classical point particles in nature, and the 
theoretical status of such objects is problematic at best, it is important to ask how any such 
prescription is constrained by physics. 

Our view is that the force law should reflect the force on an extended body coupled to a 
scalar field in the limit of small spatial extent. In particular, fix the background spacetime 
and consider a family of extended bodies and corresponding scalar fields parameterized by e, 
the spatial size of the bodies. For each body in the family, we define a center of mass world 
line z{t) (e.g., by the methods of Beiglbock [T7|] ) and calculate the charge, q, and mass, m, 
of the body with respect to this world line,[] as well as the force /|[e] exerted by the scalar 
field on the body. (For the definition of the force exerted on a small body by a field to 



1 Because the scalar charge density is a scalar quantity, the total charge that one calculates for an 
extended body depends upon the spacelike surface used to slice the body. This is in contrast to 
electromagnetism, where the charge density is the time component of a conserved vector field and 
the total charge is independent of slice. 
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which it is coupled, see Quinn and Wald 0.) We further require that m and q vanish as e 
goes to zero. For such a one-parameter family, Quinn and Wald argue that it is possible 
to specify some set of conditions on the internal structure and composition of the extended 
bodies such that the limit of /g[e] for small e is independent of their internal details. We 
would like our expression for /g to correctly the describe the order q and q 2 contributions 
to /g [e] which are independent of the internal structure of the body under these conditions. 
(Other corrections which arise from the internal structure, such as multipole effects and spin 
effects, have been derived elsewhere and should simply contribute additively at this order.) 

Unfortunately, the limit described above is quite delicate, and the task of specifying 
conditions to ensure its convergence appears to be formidable. (The analysis of Dixon [TJ[ 
demonstrates the degree of complexity which arises even without considering self-field ef- 
fects.) Nevertheless, certain properties of this limit are strongly suggested by the nature 
of the divergences in the scalar field. Following Quinn and Wald 0], we will introduce 
these properties as axioms, and then give the unique prescription for /| which satisfies these 
axioms. 

In the next subsection, we will motivate our crucial Comparison Axiom by considering the 
point particle limit described above and develop the expansions required to state the axiom. 
Then, in the following subsection, we state both axioms and give the unique prescription for 
/g that satisfies them, which is the main result of this paper. 



A. Motivation for the Comparison Axiom 

Consider a spacetime (M, g^) containing a spatially compact body characterized by 
stress-energy T^ d and scalar charge density p, a smooth Klein-Gordon field 0, and possibly 
some other set of fields which are coupled to the body, characterized by T e " fe t . The Klein- 
Gordon field satisfies the equation 

V a V a = -4vrp (1) 

with stress-energy 

T f = i-(V>VV - l -g ab g cd W c <py d <P). (2) 
Assuming that the total stress-energy is conserved, so that 

V 6 (T« 6 dy + + T e t) = 0, (3) 
then the force density exerted on the body by the scalar field is given by 

V 6 T« b dy + V b T£ t = -V b T« b = pV>. (4) 

Therefore, naively taking the point particle limit, we would expect the force on a scalar 
particle of charge q to be given by 

/s a = ?V> (5) 
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Unfortunately, this expression is meaningless as it stands, since V a diverges on the world 
line of the particle. (The situation is exactly the same with the Lorentz force law = 
qF ab Ub.) However, if we consider two points P and P along the world lines of two different 
particles in two different spacetimes (each with charge q) , and we identify the neighborhoods 
of P and P, then we might hope that, under some conditions, the difference V a — V a will 
be finite even as the two individual fields diverge. Under such conditions, it seems reasonable 
to expect that the difference between the forces on the particles will be given by the (finite) 
difference between the field gradients. That is, 



(Here, the average over a sphere of radius r, denoted by () r , is introduced to allow for the 
possibility that the r — + limit of the difference is finite, but direction-dependent.) 

Quinn and Wald give plausibility arguments which suggest that the counterpart of 
Eq. (fl) is indeed a property of the point particle limit in the electromagnetic and gravita- 
tional cases. These arguments generalize straightforwardly to the scalar case, so we will not 
give the details here. Instead, we will simply impose Eq. (H) as an axiom and investigate 
the consequences for /g. This idea will be the basis of our crucial Comparison Axiom in 
the next subsection. However, first we must find out what conditions to impose on the 
spacetimes, the world lines near P and P, and the identification of their neighborhoods in 
order to ensure that the difference in the field gradients be finite as r — * 0. In order to 
answer this question, we will now examine in detail the singularity structure of the scalar 
field in the neighborhood of the world line. 

Consider a scalar field satisfying Eq. ([I]) in a spacetime (M, g a b) with a point particle 
source 



In contrast to the electromagnetic case, the Klein-Gordon equation does not require conser- 
vation of charge. For simplicity, we shall assume throughout our analysis that the charge q 
is constant along the world line. We wish to expand <fi in r, the spatial distance from the 
world line z{r). We are primarily interested in the divergent contributions to 0, character- 
ized by the negative powers of r in the expansion, since these divergent contributions will 
determine the conditions required for convergence of the limit in Eq. (^j). It follows from 
the general theory of propagation of singularities (see theorem 26.1.1 of Hormander |19| ) 
that every solution of Eq. (|I|) which is smooth away from the world line will have the same 
singularity structure near the world line, so we are free to choose any convenient solution 
for our expansion. Later, when we wish to produce an explicit expression for /g, we will 
want to write (p m terms of the advanced and retarded solutions. Therefore, these are the 
solutions which we will analyze in the following expansion. 

Given any point x in a spacetime (M,g a b), there exists a convex normal neighborhood 
C(x) containing x [i.e. a neighborhood C(x) such that there exists a unique geodesic con- 
necting any two points within C{x)\. For x' G C(x), the Hadamard elementary solution of 
Eq. ([!]) can be written in the form @] 



/ s a -/ s a = limg(V a 0-V a 0) r . 



(6) 




(7) 
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G (1 \x,x') 



7T 



U(x, x' 



+ V(x, x 1 ) In \a(x, x') | + W(x, x') 



a(x, x') 

with corresponding advanced (+) and retarded (-) Green's functions 

G±(x, x') = 9±(x, x') U(x, x')5( a(x, x ) ) — V(x, x')9[ —<j{x, x') 



(8) 



(9) 



Here, a(x, x') is the biscalar of squared geodesic distance?] and U, V, and W are all smooth 
biscalar fields. (For an explanation of the bitensor formalism, see Dewitt and Brehme Q.) 
The scalar function 6±(x, x') is unity when x' is in the causal future/past of x and vanishes 
otherwise. 

For x near the world line z(t), let be the proper time along the world line which 
is simultaneous with x in the sense that the spatial surface £ generated by geodesies per- 
pendicular to u a at z(te) intersects x. In particular, let x lie a proper distance r along 
the geodesic generated by unit spatial vector f a at z(ts), and let z{r + ) and z{r_) be the 
intersection of the world line with the future and past light cones of x, respectively. We 
require that x be close enough to the world line that z(t^), z(t + ), and z(t_) all lie within 
the neighborhood C(x), and we denote the future and past intersections of the world line 
with the boundary of C(x) by z(T + ) and z(TL), respectively. This is illustrated in Fig. |l|. 
For the retarded field </>_, we then have 



(x) = I G-(x, x')p(x')y/—g d 4 x' 



G-(x, x') [ I q5 A (x\ z(t)) dr ) \f—g d A x' 



q G-(x, z(t)) dr 



q I + 9-[x,z(t)] U{x,z{T))8(a{x,z{T)) S j - V(x, z(r))e^-a(x, z(t))^ 



dr 



+ q G-(x,z{r))dr 

J —oo 

q n\U5{a) - V6{-a)\ dr + q [ G^dr 

JT- J-oo 



(10) 



In the last line and hereafter, we suppress the spacetime dependence for all biscalars, since 
each depends upon x in its first argument and z(t) in its second argument. For a bitensor 
A, we introduce the notation 



2 The biscalar of squared geodesic distance a(x,x') is equal to half of the squared length of the 
geodesic connecting x and x': negative for timelike separated events, positive for spacelike separated 
events, and zero for null separated events. It is defined only when there is a unique geodesic 
connecting x and x' . 
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A = ^-A{x, z(r)) = u a 'Va>A(x, z(r)). 

(IT 



We have 



dr 



do 



do\ 1 



dr = -—do = — do = & do 



drj 



(11) 



(12) 



which gives us 



{a- l U} r=T _- / Vdr 



q / G-dr. 



(13) 



We now wish to produce the corresponding expression for V a 0_. Note that the right 
side of Eq. (|13|) depends upon x in two ways: explicitly through the first argument of each 
biscalar and implicitly through r_. We have 



Vdr 



+ q V a G^dr 



{-&- 2 V a oU + o- l V a U} T=T _ + {-o- 2 oU + cT 1 ?/} V a r_ 



V a Vdr-{V} T=T _V a i 



q / V a G_ dr 



Since o(x, z{rJ)) = 0, we have 

V„{o-} T=T _ = {V a a} r=r _ 
so that 

V a r_ = {-^V^}^.. 

Therefore, we have 



{cr} T=T _V a r_ = 0, 



(14) 



(15) 



(16) 



V a <; 



Q 



{-o- 2 V a oU + o- l W a U + &- 3 oUV a o - o- 2 UV a o + <rVV a <r} 



VaVdT 



+ q VaG.dr 



(17) 



In Eqs. (|T^) and (|T7D, we would like to combine the integrals which appear on the right 
side. For T_ < r < r_, we have G(x, -z(r)) = — V(x, z(t)). Furthermore, since V is a smooth 
biscalar, 



V dr = lim 



(18) 



and 
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V a VdT = lim / VaVdr. (19) 
e— o y T _ 

Therefore, combining the integrals, we have 

0_ = q {a^U}^ + \imq f GL dr. (20) 

i/ — oo 

and 



+ lim q I V a G_ dr 



(21) 



In order to investigate the singularity structure of 0_ and V a 0_ near the world line, 
we need expansions for the expressions in brackets on the right sides of Eqs. (|20]) and fl2T|) 
which are valid to 0[r°]. (The integrals in these equations make smooth contributions to 
the fields.) The required small distance expansions for U, V, a, and their derivatives can 
all be found in DeWitt and Brehme j| or derived straightforwardly from expressions given 
therein. Switching the roles of the primed and unprimed indices for notational simplicity 
and including the corresponding results for the advanced field, </>+, we have 



x') = q (V 1 - -a a r a \ ± limg / G±(x', z{t)) dr + 0[r) (22) 
V 2 J e^O J , 



and 



V a 4±{x') = qg a >a (-r- 2 r a - -r^V + K~\a b r b )r a - ^-{a b r b ) 2 r a + ^(a b f b )a° 

- ^R bdce u b u c f d r e r a - ^a 2 f a - ^R bc f b r c r a + ^(d b f b )u a 

+ ±^R bc u b u c r a + ^R bc u b f c u a + ^R a cbd u b u c f d ± ~aV =F U a 

=F -R bc u b u c u a + -R ab r b - —Rr a =f -R ab u b ± — Ru a 
6 6 12 6 12 

/■rhoo 

±limg / V a G ± {x',z(r))dT + 0[r}, (23) 

where g a > a is the bivector of geodetic parallel displacement, defined by DeWitt and 
Brehme [|J. 

We began this calculation in order to investigate what conditions we need to impose on 
the spacetime neighborhoods and trajectories of scalar particles in different spacetimes and 
on our identification of these neighborhoods in order to ensure that the subtraction of field 
gradients in Eq. (§) is finite, and Eq. (|23|) provides the answer to this question. Since the 
divergent terms in Eq. (f23|) depend only upon the four-velocity and four-acceleration of the 
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particle (and not, for example, on higher derivatives of the motion or the local curvature), 
the subtraction in Eq. @ will be finite as long as the magnitudes of the four-accelerations 
of the two particles are equal and we identify the local spacetime neighborhoods in such a 
way that the four-velocities and four-accelerations, the geodesic distances from the world 
lines, and the parallel transport defined by g aa ' all coincide up to 0[r°]. Given points P 
and P on two world lines such that a a a a = a a a a , we can achieve this by identifying the 
spacetime neighborhoods of P and P with their respective tangent spaces Tp and Tp via 
the exponential mapQ, and then identifying Tp and Tp via any linear map which takes u a to 
u a and a a to a a . Under this identification, it is clear that four- velocities, four-accelerations, 
and geodesic distances will coincide exactly, so we need only check that parallel transport 
will also agree up to the appropriate order. 

One way to see this is to write out Eq. (|23| ) explicitly in coordinates adapted to our 
identification map, so that each point in the neighborhood of P is mapped to the point 
with the same coordinates in the neighborhood of P. (Using such coordinates, our map 
identifies a vector field in the neighborhood of P with the vector field in the neighborhood 
of P having the same coordinate components.) One such coordinate system is Riemann 
normal coordinates^. In these coordinates, the coordinate components of g a ' a are given by 

g*p = g a( 3 + \r 2 R aiP5 Pr s + O [r 3 ] . (24) 

(We have dropped the primed indices completely since expression relates components rather 
than tensors.) Comparing this to Eq. ([23D, we see that g a p simply acts as the identity at this 
order in r. [The term — (l/6)qR a y/3sr' 1 r s r a , which arises from the multiplication of the r~ 2 
term in Eq. ( p3| ) and the r 2 term in Eq. ([24]), vanishes by the symmetries of the Riemann 
tensor.] Therefore, the divergent terms will indeed cancel under the identification we have 
described. This provides the basis of our crucial Comparison Axiom in the next subsection. 

B. The axiomatic approach 

We are now prepared to give our prescription for /g , the total scalar force acting on the 
particle. We have seen that the subtraction of field gradients in Eq. flf^) will be finite as 
long as the two particles' four-accelerations have the same magnitude and we identify the 
spacetime neighborhoods via the exponential map as described above. We now elevate this 
property to the status of an axiom that any prescription for /g must satisfy. 
Axiom 1 (Comparison Axiom) Consider two points, P and P, each lying on time- 
like world lines in possibly different spacetimes which contain Klein-Gordon fields and <j) 



3 The exponential map identifies v a £ Tp with the spacetime point which lies unit affine parameter 
along the geodesic generated by v a . 

4 In order to construct Riemann normal coordinates for a neighborhood of point P, identify points 
in the neighborhood with points in Tp via the exponential map, and then pick any orthonormal 
basis for Tp. 
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sourced by particles of charge q on the world lines. If the four-accelerations of the world 
lines at P and P have the same magnitude, and if we identify the neighborhoods of P and P 
via the exponential map such that the four-velocities and four-accelerations are identified, 
then the difference between the scalar forces /g and /g is given by the limit as r — > of the 
field gradients, averaged over a sphere at geodesic distance r from the world line at P. 



Since the Comparison Axiom requires only that the four-accelerations of the particles 
agree, we now need only fix the dependence of /g on acceleration in some arbitrary spacetime 
in order to uniquely determine /g. Motivated by the time- reflection symmetry of the half- 
advanced, half-retarded solution for a uniformly accelerating trajectory in flat spacetime, 
we impose the following axiom, which should be familiar from electromagnetism. 
Axiom 2 (Flat spacetime axiom) If (M, g a b) is Minkowski spacetime, the world line is 
uniformly accelerating, and is the half-advanced, half-retarded solution, (f> = \{4>+ + 4>-)i 
then /g = at every point on the world line. 

We will now show that, if there exists a prescription for /g satisfying these two axioms, 
it must be unique. Consider a point P on the world line of a scalar particle of charge q in 
some spacetime, and let the particle have acceleration a a at point P. Let /g and g$ be two 
prescriptions for the scalar force, both satisfying the axioms given above. Now consider a 
uniformly accelerating particle with the same charge q and the same acceleration a a in a flat 
spacetime (M 4 , r] a b), and construct the half-advanced, half-retarded solution = \{<p+ + 0-) 
for this particle. By our second axiom, we know that /g = g§ = at every point P along the 
world line of this uniformly accelerating particle. Therefore, identifying the neighborhoods 
of P and P as in the Comparison Axiom above, we have 



/s a - 9s = (/s a - is) - (9s ~ ~9s) = limg(V a - V a 0) r - lima(V a - V a 0) r = 0. (26) 

r— *0 



This argument establishes uniqueness, but it also demonstrates existence by providing 
a prescription which is guaranteed to satisfy the axioms. Namely, given a point P along 
the world line of a scalar particle with charge q in any spacetime, we simply construct the 
half-advanced, half-retarded solution for a uniformly accelerating particle in flat spacetime 
with the same charge and acceleration. The scalar force /g is then given by 



This is the prescription for the total scalar force which we set out to find at the beginning 
of this section. 

Writing as = 0i n + </>_, we can use Eq. (p3|) to turn this prescription into an explicit 
formula for /g. The result is 



/ s a -/ s a = limo(V a 0-V a 0) r 



(25) 



/s a = limo(V a 0-V a 0) r . 



(27) 




1 

~ 12 



a 



) 




(28) 
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This expression, which is the main result of the paper, allows us to calculate /g for any 
trajectory z(r) in any spacetime. As stated at the beginning of the section, the physi- 
cal significance of this expression is that it should correctly describe the order q and q 2 
contributions to the force on a nearly spherical extended body in the point particle limit. 



III. THE EQUATIONS OF MOTION 

We now wish to consider the special case in which no non-scalar forces are present, so that 
the evolution of the world line z(t) is determined by the scalar field. In the next subsection, 
we derive equations of motion for z(r) in this case. Then, in the following subsection, we 
explore one of the consequences of these equations of motion: that the mass of particle varies 
with time. 



A. Derivation of the equations of motion 

Consider once again the extended body described in Sec. [IT]. In the absence of non-scalar 
fields, conservation of stress-energy dictates that 

V 6 T b a 6 dy = -V b Tf. (29) 

According to the arguments of Quinn and Wald J7J, in the point particle limit, the center of 
mass world line z{r) will therefore satisfy 

(J 771 

u b V b (mu a ) = —u a + ma a = / s a , (30) 

where /g is the limiting force we derived in Sec. [TI[ Inserting our expression for /g from 
Eq. (p8|) and separating the components parallel to u a and perpendicular to u a , we have 

a a = -(fg + u a g bc u b fi) 

777 

- -(V a in + wVV fe in ) + - (l(a a - a 2 u a ) + \{R ah u b + R bc u b u c u a ] 



\ T HI 1 V T 111 /' lO\ / /"» 

m m \ 3 o 

a 2 , 



+ lim^- J {V a G^ + u a g bc u b V c GJ)dT' (31) 
and 

dm 1 f T ~ e 

-T = ~fgu a = -qu a V a ct> in - -q 2 R - lim q 2 / u a V a G^ dr'. (32) 
dr 12 e^o 

We now note three important features of these equations. First, for each point along the 
world line, the integrals in these expressions represent that portion of V a 0_ which arises 
from source contributions interior to the past light cone of the point. This contribution to the 
force, often called the "tail term," is a direct consequence of the failure of Huygen's principle 
in curved spacetime, and can be understood as the result of scalar radiation backscattering 
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from the background curvature and re-intersecting the particle world line. The presence 
of this tail term is the primary obstacle to applying these equations in physically realistic 
situations, since most methods for calculating the retarded field of an arbitrary world line 
irretrievably mix the tail and non-tail portions of the field. 

Secondly, we can provide further insight into the nature of Eq. (|32| ) by tracing the origin 
of the Ricci scalar term in the Hadamard expansion of the field given in Sec. |I[ This term 
arises directly from the {V} T=T _V a r_ term in Eq. (]I4]). In particular, we have 

IimG_(«(r),z(r / )) = ii2, (33) 



so that we can rewrite Eq. fl32|) as 

dm 



-<KV a (0 in + tai i), (34) 



where ta ii is defined by 



'tail 



■»0 



limg / G^dr'. (35) 



The implications of Eq. ( [34]) for global energy conservation are explored by Quinn and 
Wald 0. 

Finally, owing to the presence of the Abraham-Lorentz a a term, these equations share the 
unphysical "runaway" solutions which have been discussed thoroughly in the electromag- 



netic case. (See Jackson [[Uj for one such discussion.) In order to interpret these solutions, 
it is important to remember that we view the force law given by Eq. fl28|) as an approximate 
expression for the force on an extended body, valid to C% 2 ], rather than a fundamental de- 
scription of a point particle. Therefore, we can eliminate these unphysical solutions through 
the reduction of order technique. This technique is discussed in detail by Flanagan and 
Wald |2(| , but the basic idea is simple. Recall that we wish Eq. ([H]) to describe the limiting 
motion of a one-parameter family of extended bodies in which both the charge and the mass 
of the bodies vanish as the parameter goes to zero. For concreteness, let us assume that the 
charge and mass are given by q = ae and m = be. In order to apply the reduction of order 
technique to Eq. (j3~ID , we simply insert the entire right side of the equation in place of a a in 
the a a and a 2 u a terms and discard any resulting terms which are 0[e 2 } or higher. The result 
is 



^(V^in + M^VVin) 

m 



6 m \ in x ; m z J 

1 2 2 pr—e 

+ --(R ab u b + R bc u b u c u a ) + lim — / (V a G_ + u a g bc V b G.u c ) dr', (36) 
6 m <^o m J _ oc 

which is free of the unphysical runaway solutions. 
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B. Time variation of the mass 



In stark contrast to the electromagnetic case, /g includes contributions which point 
along the four- velocity of the particle, resulting in a time-varying mass. This is not a special 
feature of the self force, nor of curved spacetime. Rather, it reflects a fundamental difference 
between the two continuum theories. Consider a small body in Minkowski spacetime with 
a center of mass world line z(r). The rest mass of such a body is given by 



m 



u aT hody ebcdei (37) 

s 



where u a is the four- velocity of z(t) (defined away from the world line by global parallelism), 
£ is the surface perpendicular to u a , and e a b c d is the volume element compatible with the 
(flat) metric. Therefore, we have 



dm d 



dr~~ dT.k ahody 



a bodv We 



/ £w[ u a ^"body e bca!e] 

- [ u a V b Tg b ody w c e cdef , (38) 



where w a is the vector field which connects successive time slices E(t). For a body coupled 

pab 
" body 



to a scalar field, we have V&T^ d = —V b T^ b = pV a 0, so that 



dm f 

= - j pu a V a (j)w c e cdef , (39) 

which is clearly, in general, nonvanishing. By contrast, in the electromagnetic case, we have 
VtJ^dy = -V 6 T^ = F ab Jb , so that 

^ = - j^u a F ab 3b( pw c e cdef . (40) 

For typical models of charged matter, j a and u a will become collinear as we take the point 
particle limit, and dm/dt will vanish. 

Perhaps because it is tempting to generalize from the more familiar electromagnetic case, 
this time variation of the mass in the scalar case has largely been ignored in the literature. 
Some authors use the equation of motion ma a = qV a 4> (e.g., Shapiro and Teukolsky [fjj). 
This equation is clearly inconsistent, and therefore in general has no solutions, since a a is 
perpendicular to the four-velocity while V a 0, in general, is not. Others explicitly project 



V a perpendicular to the four-velocity as in Eq. (|3T|) above in order to obtain the acceleration 
of the particle, but then simply ignore the component of V a which points along u a and 
assume that the mass is constant (e.g., Ori flOp. While such an equation of motion is 
mathematically consistent, it violates global conservation of stress-energy. (See Quinn and 
Wald 0.) 
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In the discussion above, we have motivated our point particle equations of motion by 
imposing local stress-energy conservation on continuum matter and taking the point particle 
limit, using our axioms to extract the appropriate finite part of the divergent fields. The time 
variation of the mass arises as a direct consequence of this local stress-energy conservation. 
In the literature on point particles, one sometimes sees an alternative derivation which makes 
no reference to the continuum theory. Instead, the author defines an action for the point 
particle system and then formally minimizes this action with respect to variations of the 
fields and the world line in order to produce equations of motion. For completeness, we give 
such a derivation here, paying particular attention to the time dependence of the particle's 
mass. 

Fix a globally hyperbolic spacetime (M, g ab ) and two Cauchy surfaces for the spacetime, 
C\ and C2. Let be a smooth scalar field and z{j) be a smooth world line in the region V 
between C\ and C 2 . We fix the value of and the position of z(r) on C\ and C 2 and define 
the action, S, as 

S = J (^V a 0V fe 0) + l - J mg ah u a u b 5\x - z{r)) dr + j q<j>5\x - z{r)) dr\ e abcd . (41) 
Formally minimizing this action with respect to variations of 0, we arrive at 

V a V a = -4vr / q5\x - z(r)) dr, (42) 



while minimization with respect to variations of z{r) yields 

dm 

—u a + ma a = qV a (f). (43) 
dr 

These are the same equations we arrived at by considering the point particle limit of the 
continuum theory. Of course, here we have assumed and z{r) to be smooth in order to 
define the action, while the solutions of Eq fl4"2"|) are clearly distributional. Therefore, no 
solutions of these equations exist. However, we may view this as a formal derivation of our 
equations from an action principle. 

Note that, if we had assumed from the outset that m was constant, the only change to 
the equations would have been to set dm/dr = in Eq. fl4"3"|). Clearly, the resulting equation 
is inconsistent, since V a does not, in general, point along the four- velocity. Still, one might 
wonder, despite the stress-energy conservation arguments given above, whether the above 
action can be modified to produce the equation of motion 

(1 777 

— u a + ma a = q(yy + u a g bc u b V c (j)), (44) 
dr 

since this equation would have the immediate consequence that dm/dr = 0, as in the 
electromagnetic case. Wiseman [^] has considered a large class of possible coupling terms 
and has found that, within this class, one cannot produce Eq. (|44T) without introducing a 
nonlinear coupling on the right side of Eq. (f4"2"D. Based on this work and the stress-energy 
considerations discussed above, we conjecture that there exists no action which produces 
Eq. (03) while preserving Eq. (B2|). 
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FIGURES 




FIG. 1. The neighborhood containing x and x' . 
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